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ABSTRACT 



In this paper, we derive the period-himinosity (P-L) relation for Large Magellanic 
Cloud (LMC) Cepheids based on mid-infrared AKARI observations. AKARI's IRC 
sources were matched to the OGLE-III LMC Cepheid catalog. Together with the 
available / band light curves from the OGLE-III catalog, potential false matches were 
removed from the sample. This procedure excluded most of the sources in the S7 
and 5'11 bands: hence only the P-L relation in the A^3 band was derived in this 
paper. Random-phase corrections were included in deriving the P-L relation for the 
single epoch AKARI data, even though the derived P-L relation is consistent with 
the P-L relation without random-phase correction, though there is a 7 per-cent 
improvement in the dispersion of the P-L relation. The final adopted N3 band P-L 
relation is 7V3 = -3.246 log(P) + 15.844, with a dispersion of 0.149. 

Key words: Cepheids — distance scale. 



X 



1 INTRODUCTION 

The mid-infrared Cepheid period-luminosity (P-L, also 
known as the Leavitt Law) relation is becoming more im- 
portant in the f uture extra-salactic di stance scale studies 
(for example, see lFreedman et al.ll2008al V This is mainly be- 
cause the effect of extinction is neRliRible in t he mid-infrared 
iFreedman et al.ll2008bl : iNgeow et all 120091 '). Other advan- 
tages of using the mid-infrared P-L relation include (see also 
the discussion in iFreedman et al] l2008bl: iNgeow fc Kanbui 
20081 : iMadore et al.ll2009l : lNgeow et al.ll2009l : iMarengo et al. 
2010|): (a) the dispersion of the infrared P-L relation is 



smaller compared to the BVI P-L relations; (b) the mid- 
infrared light curves are expected to have smaller amplitudes 
than the optical counterparts, therefore a smaller number 
(or even a single epoch) of observations are adequate to 
derive accurate mean magnitudes and P-L relations, and 
(c), metallicity may not affect the luminosity of Cepheids at 
these wavelengths (however, metallicit y may affect the ob - 
served magnitudes due to mass loss, see iNeilson et alll2009f ). 
The SAGE project (Surveying the Agents of a Galaxy's 
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Evolution. lMeixner et al.]|2006l ) has surveyed the Large Mag- 
ellanic Cloud (LMC) using the Spitzer's IRAC instrument 
at two epochs. The first data released from the SAGE 
team only include d Epo ch 1 data. This has b e en used by 
i Ngeow fc Kanbuii (|2008l ') and IFreedman et al.l (|2008bl ') to 
derive the LMC P-L relation in IRAC bands. The sec- 
ond release of the SAGE data contains both the Epoch 1 
and Epoch 2 data, and were used to i mpr ove the IRAC 
band P-L relations in lNgeow et al. I (|2009l l and lMadore' et al.l 
{2009), respectively. The main difference between the Freed- 
man/Madore group and our group is that the y matched the 
SAGE catalogs to the LMC Cepheids from iPersson et al.l 
HqoI), while we adopted the Cepheid catalogs released 
from t he OGLE team (Opt i cal Gravitationa l Lcnsing Exper- 
iment, [UdalskreraO [1991; [Sosiyisk^^ The two 
approaches lead to a discrepancy found in the slopes of the 
IRAChand P-L relations, in a sense that the slopes found by 
Freedman/Madore group are steeper than the slopes found 
by our group (see Table [1] for a comparison of the slopes 
found by these two groups). Nevertheless, the SAGE cat- 
alogs used by both groups do not contain information re- 
garding the time of observation. Hence, random-phase cor- 
rections cannot be applied to the single or double epoch 
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Table 1. Comparison of the P-L slopes in Spitzer's /i?j4C bands. 



Band Ngeow et al. (2009) Madore et al. (2009*) 



3.6/jm 


-3.253 ±0.010 


-3.40 ±0.02 


4.5/im 


-3.214 ±0.010 


-3.35 ± 0.02 


5.8/^m 


-3.182 ± 0.020 


-3.44 ±0.03 


S.O^jm 


-3.197 ±0.036 


-3.49 ±0.03 



SAGE data, in order to derive the mean magnitudes in IRAC 
bands. 

Independent of the S AGE project, the AKARI satel- 
Ute (|Murakami et al]l2007l ) has s urveyed the LMC u sing the 
on board IRC (Infrared Camera, lOnaka ot al."2007') instru- 
ment. Initial results were published in Ita ct al. (200^). The 
wavelength coverage of AKARPs IRC instrument is similar 
to Spitzer's IRAC bands. Therefore, the main purpose of 
this paper is to construct P-L relations based on AKARI 
observations, and attempt to resolve the discrepancy in the 
slopes of the IRAC band P-L relations calculated by our 
group and the Madore/Freedman group. The j4 A'.Ai?/ inter- 
nal catalog includes the time of observation: this allows us 
to obtain the phase information of the matched Cepheids: 
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Figure 1. Distributions of the separation between the matched 
AKARI sources and the input OGLE-III LMC Cepheids. 



mod( 



t-to. 



(1) 



where t is the time of observation, to is the time of maxi- 
mum light (either in F or 7 band) and P is the pulsational 
period. Hence £ [0, 1] represents one pulsation cycle of the 
Cepheid. The phase information allow us to apply random- 
phase corrections and derive mean magnitudes: such an ap- 
proach has not been applied to the SAGE data. In Sec- 
tions 2 and 3 we describe the data used in this study and 
our adopted methodology for random-phase correction. We 
present the results and discussion in Section 4, and followed 
by the conclusion in Section 5. Extinction correction is ig- 
nored in this paper since it is negligible in the mid-infrared 
bands. 



2 THE DATA 

The AKARI cutnlog provides photometric data at 3 (A'^S), 
7 (57), 11 (Sll), 15 (il5), and 24 (L24) microns together 
with their time of observation. Photometry i s on the IRC- 
Vega magnitude system, which is defined in iTanabe et al.l 
l|2008l ). Details of the data red uction and catal og compi- 
lation processes are described in llta et al.l (|2008l ). and will 
not be repeated here. The coordinates given in the AKARI 
catalog are calculated by matching detected point sources 
with corresponding 2MASS sources. If matching with the 
2MASS catalog is unsuccessful (such cases usually occur in 
L15 and L24 images), then w e use the SAGE point source 
catalog (|Meixner et al.l |2006| ) as the positional reference. 
The root-mean-squares of the residuals between the input 
2MASS/SAGE catalog coordinates and the fitted coordi- 
nates are smaller than 1.2 arc-second for N3, 2.6 arc-second 
for S7 and 5*11, and 2.9 arc-second for L15 and Z/24, respec- 
tively. The AKARI catalog coordinates should be accurate 
to that extent relative to the 2MASS and SAGE catalog 
coordinates. 
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Figure 2. The P-L relations constructed from all of the matched 
AKARI sources to the input OGLE-III LMC Cepheids, with a 
crude division of the matched sources that has separation greater 
(open circles) and smaller (filled triangles) than 1 arc-second. Er- 
ror bars are omitted for clarity. 
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We matched t he AK ARI's IRC point source catalog 
given in llta et al.l l|2008l fl to the OGLE-I I I fund amental 
mode Cepheid catalog from ISoszvriski eraP (|2008l ). Due to 
the smaller area co verage of the AKARI survey (see Figure 
1 of llta et al.ll20oj ) as compared to the SAGE project, the 
number of matched sources is reduced to 537, 226 and 83 in 
the AKARI's N3, S7 and 5*11 band, respectively. Figure [1] 
shows the distribution of the separation between matched 
AKARI sources and the OGLE-III LMC Cepheids. The 
numbers of matched sources with separation greater than 
1 arc-second are 100 (18.6%), 124 (54.9%) and 55 (63.3%) 
in the A'^3, 57 and 5*11 band, respectively. The correspond- 
ing P-L relations from all the matched sources are presented 
in Figure (2] This Figure shows that there is a well defined 
sequence of P-L relation in 7V3 band. The P-L relation se- 
quence is not obvious in the 5*7 band and almost disappears 
in the 511 band. A fraction of the matched sources devi- 
ate from the expected P-L relation sequence, with most of 
them having a separation greater than 1 arc-second. They 
are clearly the false matched j4_ftrj4_R/ sources to the OGLE- 
III Cepheid catalog. 



2.1 Selection of Data 

Since the phase information for the matched j4_K'j4 7?/ sources 
can be obtained from equation (1), we can construct the 
colour of the matched sources using the f ull I band light 
curves from OGLE-III catalog given in ISoszviiski et al.l 
l|2008l) . We match the phased / band magnitude that is clos- 
est to the phase of the j4i("j4_R/ magnitudes in all three bands, 
and construct the I{(j)) — AKARI{<I)) colour at a given phase 
<f) for a particular matched source. The phase difference be- 
tween the I band magnitudes and the AKARI sources is 
always less than 0.03, with majority of them smaller than 
0.003. We then plotted the /((/>) - AKARI{(I)) colour as a 
function of phases and separations in Figure [3] 

It can be seen from upper left panel of Figure |3] that 
colours with 7(0) — N3{<l)) < 2.0 show a periodic behavior. 
This is because most of these AKARI sources are indeed 
matched to the Cepheids, and hence this periodic behav- 
ior originates from the pulsational behavior of Cepheids. 
In contrast, the colours with IJA) — S7(0) < 2.0 did not 
show a clear periodic behavioiQ, and the periodic behav- 
ior is totally absent for 1(0) — 511(0) < 2.0 colours. The 
lower panels of Figure [3] suggest that the matched AKARI 
sources with /((;/>) — AKARI{<I)) > 2.0 in general have a 
larger separation: this further supports the conclusion that 
these AKARI sources are falsely matched to the OGLE-III 
Cepheids. Therefore, we only consider the matched AKARI 
sources in A'^3 band in this paper with the following selection 
criteria (based on Figure |3I : 



Separation < 0.7arc-second, 



^ A revision of this catalog is currently in progress (Kato et al., 
2010, in preparation), and will be published in the near future. 
Preliminary analysis shows that the difference of the band 
photometry is negligible between the current and revised catalogs. 
^ Even if the I(<f>)—S7{4>) < 2.0 colours show a very weak periodic 
behavior, the number of possible truly matched AKARI sources 
is small for a meaningful analysis of the P-L relation in S7 band. 



0.5 < I{(l>) ~ N3{(f)) < 1.7. 

These selection criteria are represented as a dashed box in 
the lower left panel of Figure O Figure!?] shows the random- 
phase colour magnitude diagram (CMD) for the matched 
AKARI sources. Sources that satisfy the above selection 
criteria appears to define an instability strip occupied by 
Cepheids. 



3 THE RANDOM-PHASE CORRECTION 

The methodology for random-phase correction in the 
JHK bands ha s been dev eloped by Nikolacv ot al. (2004}), 
ISoszvilski et al.l (|2005f ) and lMaiae"ss et al.l ([2008) in order to 
derive mean magnitudes from single epoch 2MASS data. All 
of these random-phase correction m ethods have the i r own 
advantages and disadvantages. The iNikolaev et all (120041 ) 
method is straight forward to apply, but their approach 
ignores th e contribution of the amplitudes in individual 



Ce pheids llSoszvnski et "aL l l2005lV_ The technique dev eloped 
bv ISoszviiski et all (|2005l ) and iMaiaess et al.l (120081 ) takes 



the amplitudes of the light curves into account. However, 
their method also relies on the existing light curve templates 
(for example in JHK bands) and/or the amplitude ratios 
between the JHK band and the V (or 7) band. 

In our case, the amplitude ratio between A'^3 and V (or 
7) band is an unknown. Furthermore, there are no Cepheids 
with full light curves in the A''3 band to develop the light 
curve template. We therefore need a modified approach to 
derive the random phase correction for our single-epoch A'^3 
data. 



3.1 The Methodology 

Following ISoszvhski et al.l (|2005l ). we define a normalized 



light curve at phase 4> as: 



m{(f))— < m > 
A ' 



(2) 



where < m > and A is the mean magnitude and amplitude 
of the light curve in a particular band, respectively, and is 
the phase calculated from equation (1). Th e function T((j>) 
can be modeled as an n-ord er Fourier series (|Nikolaev et al.l 
|2004 ISoszvnski et"al]l2005l ) : 



T(0) = 



aj cos{2t7 jtj)) + bj sm{2nj(j))]. 



(3) 



Re-arrange equation (2) such that m =< m > +A x T{(f>). 
If we denote the amplitude ratio between A'^3 band and I 
(or V) band as Am2./Ai — C, then we have m —< m > 
+AiC X T{(j}). The amplitude ratio C can be absorbed into 
the Fourier coefficients aj and bj in equation (3). Further- 
more, the mean magnitude < m > can be replaced by the 
definition of P-L relation, < m alog(P) + (3. We then 
obtain a linear equation which can be solved with a standard 
least squares method: 

m = alog(P)-|-/3 
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Figure 3. The I {(j)) — AK ARI {(p) colours of the matched AKARI sonvces to the OGLE-III Cepheids as a function of phases (top panels) 
and separations (bottom panels), for the A^3 (left pane), S7 (middle panel) and 511 (right panel). The dashed lines in lower left panel 
are the adopted cuts to the data (see text for details). Error bars are omitted for clarity. 
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Figure 4. Left panel shows the random-phase CMD for the I{<f>) — N3{<f>) colour. The dashed lines are the adopted colour cut at 
0.5 < — N3{4') < 1-7. Right panels show the random-phase CMD for the I{<j>) — S7{!j>) and I{4>) — Sll{(j>) colours. Error bars are 
omitted for clarity. 



+Ai X ^^[cij cos(27rj0) -I- b'j sm{2-K j(f>)], 



(4) 



where a'j = Cuj and b'j = Cbj . Dropping the Aj and C term 
from equation (4), we recover the e xpres sion that is similar 
to the one given in iNikolaev etahl (|2004l ): 



m = Q log(P) + P + "^[cj cos{2tv j<j)) + bj s,in{2-Kj(f>)]. (5) 



3.2 Test with the K Band Data 

The performance of equation (4) and (5), with the random- 
phase correction term, to solve for the P-L relation can be 



tested with K band data taken from iPersson et al. I (|2004h 
We first matched the 88 Cepheids (afte r removing HV 883 



HV 2447, HV 2883 and HV 12765, as injPersson et al.l 



with OGLE-HI Cepheid catalog from lSoszviiski et al 



200. 



200i 



that has both the / band amplitudes and to information. 
There are 46 Cepheids found to be matched between these 
two catalogs. Each of these Cepheids has ~ 19 to ~ 30 A" 
band data points per light curve. 

For each Cepheid, we randomly select one data point 
from the K band hght curves to represent the single epoch 
observation. We then fit the randomly selected data in three 
cases: (a) P-L relation without random-phase correction 

(i.e., -^random ^ 

log[P] -I- P only); (b) P-L relation with 
random-phase correction from equation (4) that includes an 
amplitude term; and (c) P-L relation with random-phase 
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Figure 5. An example of the comparison for the derived P-L 
relations in various cases. Upper left panel shows the P-L relation 
from using the mean magnitudes. Upper right panel shows the 
fitted P-L relation if the single epoch K band data is randomly 
selected, without any random-phase correction. Lower left and 
right panels show the resulted P-L relations if equation (5) and 
(4) is employed, respectively. In both cases an n = 2 order Fourier 
expansion is adopted. The randomly selected data presented here 
only represent one of the many runs we have tried. 
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Figure 6. Comparison of the fitted T{<t>) function for the data 
presented in Figure [S] without (upper panel) and with (lower 
panel) the amplitude term. Kjiand is the randomly selected sin- 
gle epoch data, and Kpj^pix = alog(P) + /3, where a and /3 are 
fitted results from equation (4) or (5). 



correction from equation (5) without the amplitude term. 
Note that there is no extinction correction and/or pho- 
tometric transformation applied to these K band magni- 
tudes, since the purpose here is not deriving the K band 
P-L relations, but to test the random phase correction. We 
also fit the P - L rela tion using the mean magnitudes from 
IPersson et all (|2004 ) for these 46 Cepheids in order to com- 
pare with the P-L relations derived from single epoch data. 
Figure [5] presents the comparison of the P-L relations in 
these four cases from one run of the randomly selected sin- 
gle epoch K band data. Figure [6] shows the corresponding 
fitted T{(j)) function without and with the amplitude term. 



Table 2. Summary of the fitted Gaussian to the histograms for 
P-L slopes. 

Without Amplitude With Amplitude 



n 


ma 




mo 




1 


-0.005 


0.057 


-0.007 


0.049 


2 


-0.007 


0.061 


-0.010 


0.052 


3 


-0.004 


0.064 


-0.009 


0.058 


4 


-0.004 


0.071 


-0.010 


0.059 


5 


-0.004 


0.076 


-0.011 


0.067 


6 


-0.005 


0.079 


-0.012 


0.071 



Table 3. Summary of the fitted Gaussian to the histograms for 
P-L zero-points. 





Without Amplitude 


With Amplitude 


n 


ma 


SG 


ma 


SG 


1 


0.013 


0.062 


0.015 


0.054 


2 


0.014 


0.066 


0.018 


0.059 


3 


0.011 


0.074 


0.017 


0.060 


4 


0.011 


0.080 


0.019 


0.069 


5 


0.011 


0.079 


0.021 


0.074 


6 


0.011 


0.091 


0.023 


0.079 



We have repeated the above test procedure 5000 times 
to built up the statistics for comparing the P-L from ran- 
domly selected data, either with or without the random- 
phase corrections, to the 'true' P-L relation derived from 
mean magnitudes. Histograms were constructed for these 
5000 runs when performing such comparisons, and a sim- 
ple Gaussian function in the form of Aae~^-^~"^'^^ ^"^"g was 
fitted to these histograms (where A is the difference of the 
slopes or zero-points of fitted P-L relations). If only using the 
randomly selected single epoch data, without any random- 
phase correction, then ma = 0.017 and sa = 0.093 for the 
comparison of the slopes and ma = —0.011 and sa = 0.113 
for the comparison of the zero-points of the fitted P-L rela- 
tions. 

We also ran our test procedure and comparisons for n = 
1 to n = 6 in the Fourier series when applying the random- 
phase corrections. Table[2]and|3]summarize the fitted values 
of ma and sa from the histograms with different n for the 
P-L slopes and zero-points, respectively. Several results can 
be immediately seen from these Tables: 

• Bias of the P-L slopes, ma, has decreased from 0.017 
if random-phase correction is applied to the randomly se- 
lected single epoch data. Random-phase correction without 
the amplitude term, i.e. equation (5), seems to provide a 
better result. 

• There is not much improvement for the bias in the 
P-L zero-points when using the random-phase corrections. 
Furthermore, Random-phase corrections with the amplitude 
term (equation [4]) seems to increase this bias. 

• The dispersion of the histogram, sa, has been improved 
by ~ 15% to ~ 50% when the random-phase correction is 
included in fitting the P-L relation. Random-phase correc- 
tion with the amplitude term gives a smaller dispersion than 
the correction without the amplitude term. 

Overall, Table [2] and [3] suggest that random-phase correction 
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Figure 7. Histograms for the comparison of the P-L slopes (left panel) and P-L zero-points (right panel). In each cases, the fitted P-L 
relations from randomly selected single epoch data (Slopegj and ZPflt) were compare to the P-L relations from using the mean magnitudes 
(SlopcQ and ZPq). Upper panels are for the randomly selected data points only, without any random-phase correction. Middle and lower 
panels show the histograms if random-phase corrections with n = 1 from equation (5) and (4) are applied to the data, respectively. The 
dashed curves are the fitted Gaussian function to the histograms. The fitted values of rriQ and sq from the Gaussian are given in the 
upper-left corner of each panels. 



from equation (4) with n = 1 in the Fourier series gives the 
best fit results. Figure [7] displays the histograms with n — 1 
case. The improvement of including the random-phase cor- 
rection in fitting the P-L relations can be clearly seen in 
this Figure. We therefore include the random-phase correc- 
tion with n — 1 when deriving the N3 band P-L relation in 
the next section. 

We also compared the estimated mean magnitudes from 
the random- phas e correction (m — T \cj)]) to the 'true' mean 
magnitudes from lPersson et ahl (|2004l ). A Gaussian fit to the 
histograms for the difference in mean magnitudes reveal that 
there is no bias introduced from the random-phase correc- 
tion. The error of the estimated mean magnitudes is about 
0.06mag, for both of the random-phase correction that either 
do or do not include the amplitude term. 



4 THE N3 BAND P-L RELATION 

There are 338 matched AKARI sources in the A^S band left 
after the selection criteria described in Section 2.]0. We fur- 
ther remove one Cepheid with log(P) = 0.066 to avoid con- 
tamination from overtone Cepheids. To remove any possi- 
ble outliers presented in the data, we employed an iterative 
2.5(T-clipping procedure to fit equation (4) and (5), that in- 
clude the random-phase corrections, when deriving the P-L 
relation. The resulted P-L relations (for 327 Cepheids) are: 



^ The N3 photometry of these Cepheids will 
be available at the SIMBAD CDS database 
(|littp : //simbad . u-strasbg . f r /simbad/ ) once the revised 
catalog (Kato et al. 2010) is published. 



No random-phase correction: 
A3 = -3.198(±0.046) log(P) + 15.816(±0.031), a = 0.161 

From equation (5): 
N3 = -3.189(±0.042) log(P) + 15.811(±0.029), a = 0.148 

From equation (4): 
N3 = -3.193(±0.043) log(P) + 15.813(±0.029), a = 0.149 



The derived P-L relation without the random-phase cor- 
rection agrees well to the P-L relations that include the 
random-phase corrections. This further supports the sug- 
gestion that random-phase corrections may not be needed 
for deriving the mid-infrared P-L relatio ns, due to the ex- 
pecte d small amplitudes in those bands (|Ngeow fc Kanbuil 
I2OO8I I. However, the dispersion of the P-L relation (a) de- 
creases by ~ 7% if the random-phase correction is included. 
The P-L relations derived from either using equation (4) or 
(5) are almost identical, suggesting the amplitude term may 
not need to be included in the random-phase correction. 

Since the A3 band is closely matched to the Spitzer's 
3.6/im band, the slope of the A3 P-L relation is expected to 
close to the 3.6/im P-L relations. The slopes from the latest 
determined 3 . 6Atm P-L relations are: — 3.253 (±0.010) from 
■Ngeow et all (|2009h and -3.40(±0.02) from iMadore et all 
( 20091 ). The A3 P-L relation seems to be shallower than 
the 3.6/im P-L relations. Plotting out the A^3 P-L relation 
reveals that the 6 Cepheids with log(P) > 1.2 appear to be 
(systematically) fainter than what is expected from a ridge 
line of the P-L relation (see also the upper panel in Figure[2]). 
This can cause the P-L relation to appear shallo wer. Note 
that the saturation limit at A^3 band is 7.8mag l|lta et al.l 
[2OO8), which cannot account for the appearing zero-point 
offset of the photometry for these 6 long period Cepheids. If 
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Figure 10. Comparison of tlie period dis t ribut ion for the AKARI's N3 sample we used to derive the P-L relation (upper left panel) 
and the Cepheid sample in lMadore et all l|2009l . upper right panel) to the full sample of OGLE-III Cepheids. Lower panels show the 
corresponding cumulative period distribution of these samples. 




0.5 1 1.5 0.5 1 1.5 Z 

log(P) phase (0) 



Figure 8. The final P-L relations in the 7V3 band, for using the 
single epoch random phase data only (upper panel) and with the 
inclusion of random-phase correction (lower panel) . The P-L rela- 
tions derived from equation from (4) and (5) are identical, hence 
only one of them is plotted here. Open circles are the excluded 
Cepheids (see text for details). Error bars are omitted for clarity. 

we excluded these long period Cepheids, then the fitted P-L 
relations (for 321 Cepheids) are: 

No random-phase correction: 
iV3 = -3.254(±0.051) log(P) -I- 15.850(±0.034), a = 0.160 

From equation (5): 
N3. = -3.241(±0.047) log(P) + 15.842(±0.031), a = 0.147 



Figure 9. The corresponding T{(f)) function for the Cepheids 
presented in Figure [S] Upper and lower panel shows the T{ij>) 
function without and with the amplitude term, respectively. 



From equation (4): 
iV3 = --3.246(±0.047) log(P) + 15.844(±0.031), a = 0.149 

Figure [5] presents the final A'^3 P-L relations, and Figure [5] 
shows the corresponding T{(j>) as a function of phase for the 
final adopted Cepheids. Note that the slope of the A'^3 P- 
L relation without the random-phase correctio n is identical 
to the Spitzer's 3.6^m band P-L relation from lNgeow et all 
(|2009l ). which the random-phase corrections were not ap- 
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Table 4. The effect of period cuts on the A^S P-L slope. 



log(-Pcut) 



Slope 



N < log(P) > 



0.50 


-3.203 ± 0.056 


261 


0.688 




0.55 


-3.156 ±0.065 


207 


0.731 




0.60 


-3.205 ± 0.074 


163 


0.774 




0.65 


-3.231 ± 0.085 


131 


0.810 




0.70 


-3.293 ±0.113 


91 


0.869 




0.75 


-3.399 ±0.135 


75 


0.902 


m 

QJ 


0.80 


-3.524 ±0.151 


60 


0.934 


T> 


0.85 


-3.468 ±0.209 


43 


0.978 


DEC 



plied to the SAGE data. We took the P-L relation derived 
from equation (4) as our final adopted A^S P-L relation: 



iV3 : 



-3.246 log(P) + 15.844 



(6) 



4.1 Comparison to the 3.6/im P-L Relation 

The slope of 



3.246 from the A3 P-L relation is closer to the 
3.6/im P-L slope giv en in Ngeow et al., ( 2009i . —3.253) rather 
than the slope from lMadore et al.l ( 2003 . -3.40). A possible 



explana tion for this discre pancy is th e different sample siz e 
used in lNgeow fc Kanbuj f20 08) and iMadore et al.l Hooi). 
The sample in Madore et al.l l|2009l ) consist of 71 Cepheids 
with most of them have period longer than 10 days. The 
mean period of their sample is < log(P) 1.242, in con- 
trast to OGLE-III (0.622) and the A3 sample used in this 
paper (0.641). In the upper left panel of Figure [TIJl we com- 
pare the period distribution of the final AKARI's A3 sam- 
ple that we used to derive the P-L relation to the ~ 1800 
Cepheid sample from OGLE-III catalog. It can be seen that 
both samples show a consistent period distribution. This 
is further supported from the cumulative distribution of the 
two samples, as presented in lower left panel of Figure flUl In 
contrast, the period distribution from iMadore et aU (|2009l ) 
is inconsistent with the OGLE-III sample, as shown in the 
right p anel of Figure 1 101 

iNeilson et al.l (|2010l ) found that the slope of the 3.6/.im 
P-L relation b ecomes steep e r, and comparable to the 
slope given in IMadore et all (|2009l '). when a period cut 
of log(Ppiit) = 1 .05 is applied to the SAGE data from 
iNgeow et all (12009'). In Tabled we present the slopes for the 
A3 P-L relation when various period cut (log[Pcu^]) is ap- 
plied to the A3 data, together with the number of Cepheids 
left in the sample (A) and the mean period of the sample 
(< log[P] >). From this Table, the slope of the A3 P-L re- 
lation and the number of Cepheid in the sample are almost 
same as the results given in lMadore et all ()2009il . if a period 
cut of log(P(,y^) ~ 0.75 is applied. 

Since some of the AKARI sources located at the LMC 
bar region (where the stellar density is higher), could the 
problem of crowding affects the slope of the A3 P-L rela- 
tions as seen in Table |4]? The locations of the Cepheids from 
our final adopted sample are shown in Figure [TTl where the 
dashed box represents the LMC bar region. We divided our 
sample into two sub-samples: one within the bar region and 
one outside the bar region. In Table [5] the P-L slopes de- 
rived from these two sub-samples were compared, together 
with their mean period and the number of Cepheids sepa- 
rated at log(P) — 0.75. Due to infiuence of the crowding 
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o log(P)<0,75 
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RA [deg.] 

Figure 11. Location of the Cepheids in our final adopted sample, 
separated by log(P) = 0.75. The dashed box roughly outline the 
location of the LMC bar region. 



on the photometry, especially for the short period Cepheids, 
the P-L slope for the sub-sample from the bar region is ex- 
pected to be shallower. However, this is not seen from the 
results given in Table [5] which suggest that crowding within 
the bar region does not affect the AKARI's photometry. 
Furthermore, the slopes from the two regions are consistent 
with each other. We have shifted, expanded or contracted 
the bar region as shown in Figure 1111 the results still hold 
with similar slopes presented in Table [5l 

The structure of the mid-infrared P-L relation has 
been found to be modified by the prese nce of Circumstel- 
lar Enve lopes dNeilson et al.ll2009l , [201ol 'l. like those discov- 
ered by iKervella et al.l (l2006l l. It is argued that the en- 
velopes are formed by a Cepheid wind where dust forms 
at a large distance from the star and causes an IR excess 
dNeilson fc Lester! l2008l l. This IR excess changes the struc- 
ture of the mid-infrared P-L relation, by increasing the ob- 
served mid-infrared brightness of Cepheids. The IR excess 
is found to be important in two regimes, the first is for 
Cepheids with log(P) < 0.7, and for long-period Cepheids. 
The IR excess of short-period Cepheids is significant and 
makes the slope of the P-L relation appear more shallow 
than it actually is. For long-period Cepheids, the IR excess 
causes the slope to become steeper. The observed 3.6^m P-L 
relation for log(P) > 1.0 is -3.23 (|Ngeow et all 120(391 '). and 
when the IR excess is removed the slope is predicted to be 
—3.19. For increased period cuts, the predicted P-L relations 
tend to have a shallower slope. It was also noted that the 
observed P-L relation tends to steepe n with increased pe - 
riod cut and show n that the samples in iNgeow et al.l (|2009l ) 
and IMadore et aL I lioog') are consiste nt. Thi s sug gests that 
a fraction of the Cepheids in Madore et al.l (|2009 ) sample 
coul d be affected by I R exce ss due to mass loss. 

iFreedman erahl (|2008bl ) suggested that the slopes of 
the mid-infrared P-L relation should reach an asymptotic 
value that can be predicted from the slope of the period- 
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Table 5. Comparison of TVS P-L slopes from two regions. 





In the LMC bar 


Outside the LMC bar 


P-L slope 


-3.275 ±0.079 


-3.240 ±0.056 


< log(P) > 


0.627 


0.658 


'^log(P)<0.75 


135 


111 


-'Vlog(P)>0.75 


31 


44 



radius (P-R) relation. If ignored the temperature term at 
the mid- infrared, the predicted luminosity is L cc R^, and 
hence m = —2.5 log(L) oc —5log(R). Adopti ng the observed 
P-R slope of 0.68 from iGieren et all l|l999l ). then the pre- 
dicted P-L slopes in the mid- infrared is -3.40 . This value 
is closer to th e slope s found in'Preed man et ahl ll2008bl ) and 
iMadore etal] l|2009l fl However, Neilson et al.1 l|2010l ) argue 
that the temperature term cannot be ignored, since the mid- 
infrared luminosity is at the Rayleigh- Jeans tail of the black- 
body function [L oc R^T^f^). If including the temperature 
term, then the predicted P-L slopes at the mid- infrared is 
~ —3.23 (for more details, see iNeilson et allboiol ). which is 
closer to the slopes found in the A'^3 band from this paper 
and in lNgeow et al] (|2009| ). 



5 CONCLUSION 

In this paper, we matched the AKARI's IRC catalog with 
OGLE-III LMC Cepheid catalog. This allows us to derive 
the P-L relation in mid-infrared bands. We have applied a 
conservative cut on the matched sources based on the colour 
information (at random phase) and the separation between 
the matched sources in the two catalogs. Only the matched 
sources in N3 band have sufficient number of Cepheids for 
deriving the P-L relation after applying these selection cri- 
teria. 

We have employed a modified rando m-phase correc- 
tion b ased on the methodology ad opted from lNikolaev et al.l 
l|2004 ) and ISoszvhskT et al.l l|2005h . The correction function 
was modeled with a lower order Fourier series which either 
include or witho ut the amplitude ter m. The K band light 
curve data from IPersson et al.l (|2004l ) was used to test the 
fitting procedure including the random-phase correction. It 
was found that random-phase correction with n = 1 in the 
Fourier series provides the best fit results. 

The derived N3 band P-L relation with random-phase 
correction is consistent with the P-L relation without the 
random-phase correction. However, the dispersion of the P- 
L relation was improved by ~ 7% when the random-phase 
correction is applied to the single epoch A'^3 band data. The 
slope of the A'^3 P-L relation, which is independent of the 
SAGE data, i s in good agr e emen t with the 3.6pim P-L rela- 
tion found in lNgeow et al.l (|2009| ). 



* If we accept this argument, our slopes of —3.25 from the 
A'^3 and 3.6/im P-L relations are consistent with the theoretical 
P-R slope from IBouo et al.l | |l998l . 0.653 from the non-canonical 
models for LMC metaUicity at Z = 0.008). 
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